Boundary polarization in the six-vertex model 
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Vertical-arrow fluctuations near the boundaries in the six-vertex model on the two-dimensional 
N x N square lattice with the domain wall boundary conditions are considered. The one-point 
correlation function ( "boundary polarization" ) is expressed via the partition function of the model 
on a sublattice. The partition function is represented in terms of standard objects in the theory of 
orthogonal polynomials. This representation is used to study the large N limit: the presence of the 
boundary affects the macroscopic quantities of the model even in this limit. The logarithmic terms 
obtained are compared with predictions from conformal field theory. 



PACS numbers: 05.50.+q, 05.70.Np, 02.30.1k 
I. THE MODEL 

In this paper, we shall consider the six-vertex model on 
a square lattice. Originally, this model was introduced 
as a model describing the ferroelectric properties of the 
hydrogen-bonded planar crystals fit] . The hydrogen atom 
positions are specified by attaching arrows to the lattice 
edges. In the six-vertex case, the arrows are arranged 
in such a way that there are always two arrows point- 
ing away from, and two arrows pointing into, each lattice 
vertex (the so-called "ice rule"); thus, there are six pos- 
sible states at each vertex (see, Fig. Q). The statistical 
weights, a, b and c, of the allowed states are invariant 
under the simultaneous reversal of all arrows: 



FIG. 1. The vertices of the six- vertex model and their sta- 
tistical weights. 

The partition function of the model on an N x N square 
lattice is obtained by summing over all possible arrow 
configurations {C}, 



all arrows on the left and right boundaries are pointing 
outward (see, Fig. g). 



FIG. 2. The domain wall boundary conditions. 

This model was introduced in Ref. || in connection with 
the calculation of the correlation functions for exactly 
solvable 1+1 dimensional models g. It appears that 
some problems from the theory of alternating sign ma- 
trices J7|,|| and domino tilings || may be reformulated in 
terms of this model. 

Aperiodic boundary conditions arc of interest since 
they demonstrate the influence of the boundaries and in- 
ternal defects of real physical systems on their bulk prop- 
erties. By renormalization group method, it was shown 
that the behavior of the correlation functions near the 
surfaces and defects is quite different from the bulk be- 
havior p0| |. The exactly solvable six- vertex model with 
the DWBC provides us the opportunity to study the sur- 
face phenomena beyond renormalization group scheme. 
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where rii, i%2, and 713 are the number of vertices of type a, 
b, and c in configuration C, respectively {n\ + Ti2 + n^ = 
N 2 ). 

The six-vertex model was studied for both periodic 
(PBC) P,U and fixed boundary conditions W. In this 
paper, we are concerned exclusively with the domain wall 
boundary conditions (DWBC), namely, all arrows on the 
top and bottom of the lattice are pointing inward while 



II. THE PARTITION FUNCTION AND THE 
BOUNDARY POLARIZATION 

The partition function Z/v of the model with the 
DWBC may be represented as the determinant of an 
N x N matrix IIUO]. Though there exist several com- 
binatorial representations for this determinant |13| , |l4| , it 
has so far only been calculated explicitly for some special 
cases M . Significant progress has recently been achieved 



in studying the asymptotics of Zm as N — > oo. Un- 
der some special restrictions on the values of the vertex 
weights, the bulk free energy was calculated in Ref. |lj| 
by using the Toda equation fig] , A more general re- 
sult was obtained in Ref. O] by a reformulation of the 
six-vertex model as a Hermitian matrix model to which 
saddle point integration method was applied. A represen- 
tation convenient for the large N analysis was suggested 
also in Ref. @. 

Less is known about the correlation functions of this 
model for at least two reasons. First, the calculation of 
the correlators, in general, is a more complicated problem 
than the calculation of the corresponding partition func- 
tion. Second, the lack of translation invariance caused 
by the special boundary conditions introduces additional 
difficulties. Some correlation functions for the inhomoge- 
neous model (the model with the statistical weights de- 
pending on the position of the vertex) with special choice 
of the weights were considered in Ref. |ji| . 

The boundary polarization is the one-point correlation 
function that describes the probability for the arrow on 
the fixed lattice edge on the boundary to be pointing in 
either direction. The symmetry of the model allows us to 
consider only vertical arrows. Let us denote by xn the 
probability for the vertical arrow to be pointing down. 
Note that \n for PBC is independent of the position 
of the edge and is just a spontaneous polarization of the 
system || . Here, we shall discuss xn for the edge located 
at the lower-right corner of the lattice (in Fig. g this edge 
is dotted): 



FIG. 3. xn is calculated for the arrow on the dotted edge. 

It turns out that if the arrow on the dotted edge is point- 
ing down, then due to the imposed boundary conditions, 
the allowed vertex configuration at the bottom and the 
right boundaries is determined in a unique way. Actu- 
ally, in the lower-right corner one has the vertex of type 
c; thus, the rest of the 2N — 2 vertices are of type b, 
and the DWBC are valid for the residual N — 1 x TV — 1 
sublattice. Thus, 



Xn = cb 



2N-2 



Zn-i 



■<N 



(1) 



and the problem of studying the correlation function, xn, 
is reduced to the analysis of the ratio of the partition 
functions Zn-i/Z^. 



The boundary polarization on the arbitrary edge is ex- 
pressed in terms of linear combination of partition func- 
tions on sublattices. 



III. THE CONNECTION WITH ORTHOGONAL 
POLYNOMIALS 

Due to the multiplicativity of the partition function, 
one can set c = 1 without loss of generality. Hence, the 
model is characterized by only two parameters, a and 
b. In Fig. 4 the phase diagram on the (a, b) plane for 
the model with PBC is plotted (cf., Fig. 8.5 of Ref. {§). 
It may be regarded as the phase diagram for the model 
with the DWBC, in the sense that the free energy takes a 
different analytic form in the regions divided by the solid 
and dashed lines (see Ref. |TJJ for details). Naturally, 
the ground state and low temperature behavior do not 
coincide with those for the model with PBC. 




1 a 

FIG. 4. The phase diagram in terms of the weights a and 



In this and next section, we shall consider the model 
on the solid line in Fig. |j. It is convenient to use the 
following parametrization of the vertex weights on this 
line: 



a = 1/2 - x 
b = 1/2 + x 



•1/2 < x < 1/2. 



It follows from Ref. |ll|] that the partition function of 
the model is represented as the determinant of an N x N 
matrix: 



detjv M 
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(2) 



where the matrix elements of M. are given by 

M ^ = -^+k^(x), a,k = 0,l,...,N-l, (3) 
and cf>(x) = [(| - x){\ + x)] _1 . 



Matrix M. is a Hankel matrix, that is, it has constant 
entries along the antidiagonals. Its determinant can be 
expressed in terms of objects related to the theory of 
orthogonal polynomials. Let p n (£)> n > 0, be a sequence 
of monic orthogonal polynomials with weight /?(£): 



Pn(0Prn(£,)p(0 d £, = SrnnK, n, TO > 0. (4) 



Then the determinant of the Hankel matrix with the el- 
ements 



Ti-ak 



r + V(£K, a, k = 0,1,..., TV -1, (5) 



b N _ 2 bN-i- 



is equal to @ detjvH = h^b^b^ 2 
Here fto is determined from Eq. (Q), and & n are the coef- 
ficients of the corresponding three-term recurrence rela- 
tion, 

Pn+l(0 = (a n + 0P"(0 - b n Pn-l(£), « > 1, 

with the initial conditions Pa(£) = 1 and pi(£) = £ + &o- 
The function </>(x) has the following integral represen- 
tation: 

/OO 
-CJO 

Combining this representation with Eq. (g) and Eq. (^), 
one shows that detAr A4 = detjv W, where the weight p(£) 
is equal to 



p(0 



- p-^p-iei/a 



(6) 



Then, since the coefficients 6„ and /i„ satisfy the well- 
known relation b n = h n /h n -i, n > 1, we have detjv M = 
hohi ■ ■ ■ fiN-i, and the desired representation for proba- 
bihty (ffl) is pa 
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Now let us discuss briefly the general case when a and 
b are arbitrary positive constants. The partition func- 
tion is given by formulas (0) and (H), while 4>{x) should 
be changed (see, for example, ||q, |ll|]l2| ) , the weight /?(£) 
is obtained by straightforward calculations J17J], and the 
corresponding expression for ^jv+i differs from Eq. (J7|) 
by a constant. 



IV. THE LARGE TV LIMIT 

Equation (^) reduces the problem of calculation of the 
probability xn+i to the calculation of the normalizing co- 
efficient hisr. At present there exists a powerful method 
for studying the A — > oo behavior of h^ . This method is 
based on the matrix Riemann-Hilbert conjugation prob- 
lem [2l|. The corresponding asymptotic technique was 



suggested in Ref. [g^l an d worked out for orthogonal poly- 
nomials with nonanalytic weights in Ref. 123] . Using this 
technique, we obtain the main result of the paper: 



ln/ijv = 2AlnA + 2Aln 
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(8) 



where 



tp(x, A) = (-l) iv cos[2ttx(A + 1/2)] 



and the omitted terms are of the order l/[A(lnA) 3 ]. 
Note that for x = this result coincides with the one 
obtained in Ref. ]23|| . 

Having expansion iM), it is easy to find the expansion 
for \N- 
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(9) 



From Eq. (|9|), we get all increasing terms for the partition 
function Zn, 



Z N ~ C exp (/o A 2 + fiN + f 2 In A) , 
where C > is a bounded function of A, and 



(10) 



/o = In 



tt(1/4-x 2 ) 



cos(7ra) 



o, h = T2 . 



V. RESULTS AT OTHER POINTS 

At present there are several points on the phase di- 
agram, where the determinant has been calculated in 
closed form. 

(i) The "free-fermion case", a 2 + b 2 = 1 (dotted circu- 
lar quadrant on the phase diagram). On this circle, one 
has a very simple result Zjq — 1. This result can be ob- 
tained by an appropriate limit from the inhomogeneous 
partition function |g|l^|. Therefore, one has 



lnxAf = 2(A-l)lnb. 



(11) 



(ii) The "ice point", a = b = 1. All weights are equal, 
and the partition function is just the number of allowed 
configurations {C} (Sec. I). In this case one has Ref. M 
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thus, 



. ri 16 1, 16 51 51 

In Yiv — -/V In In 1 1 - 

A 27 2 27 36 TV 72 N 2 



(12) 



For Zn we have asymptotic expansion (|l0|), where /o = 
In ^S, /i = and f 2 = -^. Note that for the model 
with PBC §, / PBC = ln^J- > Jo, that is, for large N 
the number of allowed configurations for the DWBC is 
less than for periodic ones. 

(hi) The point a = b = l/y/3. In this case J7| 



Z2N+1 - r^" 



n 



(3j-l)! 



J2N 
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3 W [(2/V-1)!] 2 ' 



J2N-1, 



and we immediately get 



at, 4 1, 27 11 
^9 + 2 ^T" 187V 



(13) 



thus, /o = ln^r, /1 = and f2 = jg- The first omitted 
term in Eq. ( jl3| ) is of the order iV -2 , and depends on the 
parity of N. 

We complete this paper with the following statement. 
One can see that in expansion (|9|), which was obtained 
for a + b = 1, there exists a logarithmic term. For arbi- 
trary a and b, we state that the logarithmic terms in the 
large N expansion of In xn will take place only under the 
following conditions: weights a and b are either on the 
solid line a + b = 1, or on the dashed lines |a — b| = 1 
of the phase diagram (Fig. ||). For the model with PBC, 
the logarithmic terms could be explained via conformal 
field theory H. Similar arguments are valid for the dis- 
cussed model. The influence of the boundary condition 
changing operators J24|] should also now be taken into 
consideration. If we move away from the solid or dashed 
lines, the logarithmic terms vanish and the expansion of 
ln%jv will be in integer powers of N. Expansions (11), 
(O) and (O) indeed contain only integer powers of N, 
thus confirming the above statement. 

Finally, we would like to mention that the six-vertex 
model with any boundary conditions can be considered 
as a model for a description of interface roughening of a 
crystal surface Pq] . An important point in these stud- 
ies is the existence of exact analytical results, which are 
known for the six-vertex model with PBC jl]-|| . Wc 
believe that our analytical results for the model with 
DWBC provide one more basis for the experiments and 
simulations in this direction. 

We thank A. H. Vartanian for bringing Ref. p3| to 
our attention. This work was partially supported by the 
RFBR Grant No. 01-01-01045. 



[1 
[2 

[3] 

[1] 



[7] 

[8] 



[9 
[10 

[11 

[12 

[13; 

[14; 
[15; 

[16 

[17 
[18 

[19 
[20 

[21 
[22 

[23 

[24 
[25 



E.H. Lieb, F.Y. Wu, in Phase Transitions and Critical 
Phenomena, edited by C. Domb and M.S. Green (Aca- 
demic Press, London, 1972), Vol. 1, p. 321. 
E.H. Lieb, Phys. Rev. 162, 162 (1967); Phys. Rev. 
Lett. 18, 1046 (1967); Phys. Rev. Lett. 19, 108 (1967); 
B. Sutherland, Phys. Rev. Lett. 19, 103 (1967). 
R.G. Baxter, Exactly Solved Models in Statistical Me- 
chanics (Academic press, San Diego, 1982). 
R.J. Baxter and A.L. Owczarek, J. Phys. A 22, 1141 
(1989); M.T. Batchelor, R.J. Baxter, M.J. O'Rourke, and 
CM. Yung, J. Phys. A 28, 2759 (1995); K. Eloranta, J. 
Stat. Phys. 96, 1091 (1999); C.B. Thorn, Phys. Rev. D 
63, 105009 (2001). 

V.E. Korepin, Commun. Math. Phys. 86, 391 (1982). 
V.E. Korepin, N.M. Bogoliubov, and A.G. Izergin, Quan- 
tum Inverse Scattering Method and Correlation Func- 
tions (Cambridge University Press, Cambridge, 1993). 
G. Kuperberg, Int. Math. Res. Notices 3, 139 (1996). 
D.M. Bressoud, Proofs and Confirmations; The Story of 
the Alternating Sign Matrix Conjecture (Cambridge Uni- 
versity Press, Cambridge, 1999). 

H. Cohn, N. Elkies, and J. Propp, Duke Math. J. 85, 
117 (1996); W. Jo ckush, J. Propp, and P. Shor, e-print 



math.CO/9801068 



A. Hanke and M. Kardar, Phys. Rev. Lett. 86, 4596 

(2001); A. Hanke, Phys. Rev. Lett. 84, 2180 (2000); 

A. Hanke, M. Krech, F. Schlesener, and S. Dietrich, Phys. 

Rev. E 60, 5163 (1999). 

A.G. Izergin, Dokl. Akad. Nauk (SSSR) 297, 331 (1987) 

[Sov. Phys. Dokl. 32, 878 (1987)]. 

A.G. Izergin, DA. Coker, and V.E. Korepin, J. Phys. A 

25, 4315 (1992). 

A. Lascoux, Seminaire Lotharingien Combin. 42 Art. 

B42p, 15pp. (1999). 

C. Krattenthaler, Seminaire Lotharingien Combin. 42 

Art. B42q, 67pp. (1999). 

V.E. Korepin and P. Zinn-Justin, J. Phys. A 33, 7053 

(2000). 

K. Sogo, J. Phys. Soc. Jpn. 62, 1887 (1993); A.G. Izergin, 

E. Karjalainen, N.A. Kitanine, J. Math. Sci. (New York) 

100, 2141 (2000). 

P. Zinn-Justin, Phys. Rev. E 62, 3411 (2000). 

N.A. Slavno v Zap. Nauchn. Sem . POMI 269, 308 (2000) 

(available at cond- mat/0005295). 



V. Korepin and P. Zinn-Justin, e- print nlin.SI/000803C ; 
J. de Gier and V. Korepin, e-print math-ph/0101036. 
N.M. Bogoliubov, A.V. Kitaev, and M.B. Zvonarev, Zap. 
Nauchn. Sem. POMI 269, 136 (2000). 
A.S. Fokas, A.R. Its, and A.V. Kitaev, Commun. Math. 
Phys. 142, 313 (1991); 147, 395 (1992). 
P. Deift, T. Kriecherbauer, K. T-R McLaughlin, S. Ve- 
nakides, and X. Zhou, Int. Math. Res. Notices 16, 759 
(1997); Commun. Pure Appl. Math. 52, 1491 (1999). 
T. Kriecherbauer and K. T-R McLaughlin, Int. Math. 
Res. Notices 6, 299 (1999). 
J.L. Cardy, Nucl. Phys. B 342, 581 (1989). 
H. van Beijeren, Phys. Rev. Lett. 38, 993 (1977); E. Car- 
Ion, G. Mazzeo, and H. van Beijeren, Phys. Rev. B 55, 
757 (1997); 



